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ABSTRACT 

We present optical and near-infrared (IR) photometry and near-IR spectroscopy of 
SN 2004am, the only optically detected supernova (SN) in M 82. These demonstrate 
that SN 2004am was a highly reddened type II-P SN similar to the low luminosity 
type II-P events such as SNe 1997D and 2005cs. We show that SN 2004am was located 
coincident with the obscured super star cluster M 82-L, and from the cluster age infer 
a progenitor mass of M©. In addition to this, we present a high spatial-resolution 
Gemini-N i^-band adaptive optics image of the site of SN 2008iz and a second transient 
of uncertain nature, both detected so far only at radio wavelengths. Using image 
subtraction techniques together with archival data from the Hubble Space Telescope, 
we are able to recover a near-IR transient source co- incident with both objects. We 
find the likely extinction towards SN 2008iz to be not more than Ay ~ 10. The nature 
of the second transient remains elusive and we regard an extremely bright microquasar 
in M 82 as the most plausible scenario. 

Key words: supernovae: general - supernovae: individual(SN 2004am; SN 2008iz) - 
galaxies: starburst - galaxies: individual (M 82) - infrared: galaxies. 
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1 INTRODUCTION 

In the nuclear regions of M 82 and other nearby star- 
burst galaxies one core-collapse supernova (CCSN) is ex- 
pected to explode every 5-10 years. However, due to 
the high dust extinction in the nuclear starburst regions 
most of the actual SN e have remained undetected (e.g., 
iMattila k Meiklell200lb . Searches working at near-infrared 
(NIR) wavelengths have discovered SNe in the nuclear and 
circumnuclear regions of a few nearby sta rburst galaxies 



and luminous infrared galaxies (LIRGs) (e.g.jKankare et al] 
201I I2OO8I: fMiluzio CaDDellarollioTol : iMattila et aDl2007l : 
Maiolino et al.1 12OO2I ). Furthermore, optical searches have 



been able to discover a few SNe in normal nearby spi- 
ral galaxies with several magnitudes of v i sual extinction 
(e.g. SNe 2002hh, 2002cv: IPozzo et al.|[2006l : IPi Paola et al.1 
I2OO2I ). 

The first claim of a detection of a su pernova at op- 
tical o r NIR wavelengths in M 82 was by iLebofskv et al.l 
(| 19861 ). who reported the discovery of SN 198 6D in 2 M^m 
observations. However, it was later shown by iGehrz et al.l 
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(| 19861 ) that the source identified as SN 1986D was already 
present in their 2.2 ^tm images obtained three years ear- 
lier and has not varied significantly in brightnes s since then. 



This s ource was designated as the feature K2 b vlDietz et al 
("1986). SN 2004am in M 82 was discovered fSin ger et a" 



2004) by the Lick Observatory SN Search (LOSS) on im- 
ages from 2004 Mar. 5.2 UT but it was already present in 
their earlier images from 2003 Nov. 21.6. Being on a bright 
spot in M 82 caused their detection software to initially miss 
the new object. On 2004 March 6.9 SN 2004am was spec- 
troscopically classified as a highly reddened type II event 
showing some spectral similarity to the type II-P SN 1995V 
(|Mattila et al.L2004 ). Radio non-detections of SN 2004am 
on 2003 Nov. 14 (8.4 a nd 15 GHz) an d on 2 004 March 9 (5 
GHz) were reported bv'Beswick et al.1 (|2004 l). SN 2004am is 
still the only SN ever discovered at optical/IR wavelengths 
in M 82. 

Being a prototypical starburst galaxy at a distance of 
3.3=b0.3 Mpc, (from t he Cepheid distance to the M 81 group; 
iFreedman et al.ll200ll ) means that there is a wealth of ground 
and space based imaging of M 82. Identification of progen- 
itors in both Hubble Space Telescope ( HST) and ground- 
based images has now become frequent (jSmartt et al.ll2003 : 
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Li et al."20Q6: 'Mattila et al."2008b"; 'Maund Smarttl l2009l : 
Eraser e t al. 2011; Van Dyk et al. 2012). Furthermore, in M 



82 th e compact star popu l ation has been studied extensively 



(e.^. iLangon et al.1 lioosl; 



Bastian et aP l2007l : ISmith et all 



l2006l : lMcCradv et al.|[2003l ). Hence SN 2004am is interesting 
not just because it is the first SN discovered in M 82 in the 
optical and NIR, but also because its physical association 
with the well studied super star cluster (SSC) M 82-L allows 
a mass estimate of the progenitor star. In one recent case, 
SN 2004dj in NGC2403, the ex plosion was coincident with 
a young com pact star cluster ( Maiz-Apell aniz et al.l l2004l : 
IWang et al.L 2005). The spectral energy distributions of the 
host cluster suggest ages and masses 13-20 Myr and 2.4- 
9.6x10^ M0. If we assume that the stellar population in the 
cluster is coeval, then this le ads to e stim ates of the progen- 
itor mass of 12-15Mr:^. IVinko et al.l (|200 9) produced an im- 
proved study of the host cluster once SN 2004dj faded using 
extensive data from the UV to NIR. They estimated a most 
likely turn-off mass of between 12 and 20 M©. SN 2004dj 
was a type II-P (Vinko et al. 2006), and the estimate of the 
progenitor mass from the turn-off age is consistent with the 
mass estimates from the direct detection of SN II- P progen- 
itors on pre-explosion images (ISmartt et al.ll20()9h . 

Over the last thirty years, there have been four de- 
tected radio transients in M 82 which have been suggested to 
be SNe. The source 41.5+59.7 (jKronberg k Sramek. lll985l ) 
was seen in February 1981 at 5 GHz with a flux density 
of 7.07=b0.24 mJy, but by April of the following year had 
faded below the detec tion limit of 1.5 mJy. Another transient 
40.59+ 55.8 (jMuxlow et al. 1994 ) was detected at one epoch 
with the Multi-Element Radio Linked Interferometer Net- 
work (MERLIN) at 5 GHz with a flux density of -1.23 mJy 
in July 1992. Neithe r of these two sources were detected by 
iFenech et al. I (|2008h to a flux limit of —20 /xJy suggesting 
them to be transient sources, as opposed to sources with 
variable nature. 

More recently, a bright (—100 mJy) radio tran- 
sient in the nuc l ear re gions of M 82 was discovered by 
iBrunthaler et al.1 (|2009h in 22 GHz images from the Very 
Large Array (VLA) obtained in March 2008. Over one 
year's time the source faded by a factor of ten, show- 
ing a spectral index of about —0.8 in the VLA (1.4- 
43 GHz) observations from April 2009 indicating an opti- 
cally thick synchrotron spectrum. Furthermore, very long 
baseline interferometry (VLBI) observations from May 2008 
and April 2009 revealed a ring-like str ucture expanding by 
—23 000 km s~^ f'Brunthal er et al]|201Ql V These observations 
definitely confirmed the SN nature of the transient des- 
ignated as SN 2008iz. Monthly monitoring of M 82 with 
the 25 m Urumqi radio telescope at 5 GHz led to the de- 
tection o f SN2008iz as a flar e on top of the M 82 radio 
emission (|Marchili et aLllioiOh . The SN peaked with a flux 
of — 160 mJy on 21st June 2008, while modelling of its light 
curve yielded an accurate explosion date of 18th February 
2008. 

In May 2009, another new radio transient was discov- 
ered in the nucle ar regions of M 82 by means of MER- 
LIN observations (|Mux1ow et al.|[2009h . thus far called the 
MERLIN transient. This source reached a peak intensity 
of 0.72 mJy at 5 GHz in early May 2009, and 1.7 mJy at 
1.6 GHz on 20th May 2009. The nature of this object is stih 
unknown due to its peculiar low and apparently steady lu- 



minosity, its longevity (still observed after nine months of its 
discovery), its position with respect to the dynamical cen- 
tre of the galaxy, and tentative evidence of proper motion 
and expansion. At the moment, the most via ble explanation 
for it s nature appears to be a micro-quasar (|Mux1ow et al.l 
I2OIOI : [Joseph et all I2OII ') displaying a high ratio of radio 
to X-ray luminosity. For instance, the source Cygnus X-3, 
which is one of the strongest Galactic micro-quasars, has a 
radio luminosit y exceeding its X-r ay luminosity by an order 
of magnitude ^ Nipoti et al]l2005l ). Thus a similar scenario 
could also be plausible in the case of the MERLIN tran- 
sient. 

In this paper we present photometric and spectroscopic 
observations of SN 2004am which despite its small distance 
has remained unstudied due to the high line-of-sight dust 
extinction and the fact it occured coincident with the nu- 
clear SSC M 82-L. We make use of this coincidence to infer 
an initial mass for the progenitor. We use deep, high spatial- 
resolution imaging from Gemini-N to search for NIR coun- 
terparts of SN 2008iz and the MERLIN transient and make 
use of these to investigate their nature and the line-of-sight 
extinctions towards these sources. Finally, we discuss the na- 
ture of the radio transients in M 82 making use of existing 
radio observations and compare these with the expectations 
from the SN rate estimates. 



SN 2004AM: 
RESULTS 



OBSERVATIONS AND 



2.1 Observations and data reduction 

SN 2004am was observed with the NIR imager and spectro- 
graph LIRIS on the Wilham Herschel Telescope (WHT). 
LIRIS spectra covering zJ bands were obtained on two 
epochs, 6 March and 25 November 2004. The observations 
included in this study are listed in Table 1. The NIR imag- 
ing data were reduced using standard IRAF routines and the 
spectroscopy using the FIGARO package as a part of the 
Starlink software. SN 2004am appears coincident with the 
SSC M 82-L, therefore, all the photometry was performed 
on subtracted image s (e.g ., see Fig. 1). The ISIS2.2 pack- 
age (Alard & Luptonl ll998l ) was used to convolve the better 
seeing image to the poorer seeing one, and to match the 
intensities and background prior to the image subtraction. 
The 25 Nov. 2004 JHKs images were used as the image sub- 
traction references for the LIRIS images, as they provided a 
smooth subtraction and a further Kg epoch on 30 January 
2005 indicated that the SN was no longer detectable in the 
25 Nov Kg frame. 

In addition, we recovered optical images containing 
the site of SN 2004am from the Hubble Space Telescope 
(HST) archive. The pipeline-reduced products from the 
HST archive were used. SN 2004am was observed with the 
Advanced Camera for Surveys (ACS) HRC instrument in 
F555W and F814W filters on 5 July 2004 (SNAP 10272; PI: 
Filippenko). The site of SN 2004am was covered also with 
the ACS/WFC instrument in F814W filter on 9 Feb. 2004. 
Unfortunately, SN 2004am was saturated in this image and 
therefore we could not make use of it in our analysis. We 
attempted to detect SN 2004am in the ACS/HRC images 
using images obtained with the same instrument and filters 
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Table 1. Optical and NIR photometry of SN 2004am. The epochs are assuming the first observation of the SN by LOSS was 2 weeks 
after the explosion (see Sect. 2.3). Between 14 and 118 days unfiltered CCD magnitudes are from Singer, Pugh & Li (2004). The "und." 
labels means the SN was undetected at these wavelengths with any significance in comparison to the coincident host cluster. 



Date (UT) Epoch V unfilt I J H K Instrument/Reference 



(2004) 


(days) 














Nov. 21.5 


14 


16.0 








Singer, 


Pugh & Li (2004) 


Nov. 23.5 


16 


16.1 








Singer, 


Pugh &; Li (2004) 


Dec. 28.4 


51 


16.3 








Singer, 


Pugh & Li (2004) 


Jan. 2.5 


56 


16.4 








Singer, 


Pugh &; Li (2004) 


Jan. 17.4 


71 


16.4 








Singer, 


Pugh & Li (2004) 


Jan. 21.4 


75 


16.4 








Singer, 


Pugh &; Li (2004) 


Jan. 30.3 


84 


16.4 








Singer, 


Pugh & Li (2004) 


Feb. 5.3 


90 


16.4 








Singer, 


Pugh & Li (2004) 


Mar. 5.2 


118 


17.0 








Singer, 


Pugh & Li (2004) 


Mar. 6.9 


120 




12.95±0.08 


12.49±0.16 


12.03±0.14 




WHT LIRIS 


June 5.9 


211 








13.06±0.12 




WHT LIRIS 


July 5.3 


241 und. 




> 20.3 








HST ACS/HRC 


Nov. 25.2 


384 




und. 


und. 


und. 




WHT LIRIS 



on 7 June 2002 as reference for the image subtraction with 
the ISIS 2. 2 package. However, no significant source was de- 
tected at the location of M 82-L in the difference images and 
we conclude that SN 2004am had already faded below the 
detection limit of HST by the 5 July 2004 epoch of observa- 
tion (see Fig. 2). 

2.2 Photometry of SN 2004am 

Aperture photometry of SN 2004am was performed in the 
subtracted LIRIS images (e.g., Fig. 1) using Gaia. The JHK 
magnitudes from 2MASS for three bright field stars covered 
in the LIRIS field of view (FOV) were used for the photo- 
metric calibration. The average of the zeropoint magnitudes 
obtained was adopted for the calibration and their standard 
deviation as the photometric uncertainty. The statistical un- 
certainty in the photometry as estimated by Gaia was in all 
cases negligible. 

The residual noise from the flux of M 82-L leaves no 
obvious point source in the subtracted HST images (see Fig. 
2). Aperture photometry was performed on the residuals 
using Gaia. For this we used a 0.25" radius aperture and 
a sky annulus between 1.5 and 2.0 x the aperture radius. 
No suitable point sources were present within the ACS/HRC 
FOV and thus we adopted an aperture correction from 0.25" 
to infinite aperture from Sirianni et al. (2005). Application 
of the appropriate Vegamag zeropoint yielded m(F814W) = 
20.3, which we adopt as an upper limit for the brightness of 
SN 2004am at this epoch. 



2.3 The type and extinction of SN 2004am 

The spectra of SN 2004am at both epochs are contaminated 
by the SSC. We subtracted the high signal-to-noise spec- 
trum of M 82-L (provided bv lLancon et al.ll2008l ) from both 
epochs. This spectrum was manually scaled and the simple 
subtraction procedure resulted in an apparently flat spec- 
trum with the SSC continuum removed in the first epoch 
(6th March 2004). However the subtraction from the sec- 
ond epoch left no identifiable SN features and we conclude 



that SN 2004am had faded below the detection limit of the 
spectral signal. 

In Figure 3, the subtracted SN 2004am spectrum is 
compared to the spectra of three type II-P S Ne: 2004et 
Maguire et al 2010 ), 1995V (jFassia et al.lll998l ) and 1997D 
Benetti et al.|[200ll V All spectra were dereddened (by the 
values quoted in the original publications) and scaled to the 
continuum level at 1.3/im of SN 2004am. SN 2004am was 
dereddened with Ay = 3.7 and R = 2.4 as found for M 82-L 
by Lancon et al. (2008). SN 2004am has a similar spectrum 
to SNe II-P at the end of the plateau phase. The ejecta 
velocity from the P-Cygni troughs in the metal lines (e.g. 
Srii, Fell, and Ci) indicates that SN 2004am appears more 
similar to the low velocity and low luminosity SN 1997D, 
than the other two. Cross-correlation and fitting Gaussians 
to the centroids of the sharpest P-Cygni profiles show that 
the velocity difference between SNe 2004am and 1997D is 
negligible, while the line troughs of SNe 1995V and 2004et 
are blueshifted by 1000 it 500km s~^. Furthermore the spec- 
tral features of SNe 1995V and 2004et are visually broader 
(see Fig. 3). 

The classification of SN 2004am as a type II-P is sup- 
ported by its unfiltered light curve shown in Fig. 3, com- 
pared with the i^-band light curve s of two different II- 
P events: SNe 1999e m and 2005cs ( Elmhamdi et al.l l2003l : 
IPastorello et al.ll2009f ). The plateau points of SN 2004am are 
unfiltered CCD magnitudes from Singer, Pugh & Li (2004), 
and the upper limit at 241 days has been estimated from the 
F814W HST images. The light curves of the two supernovae 
have been scaled to the same distance as SN 2004am, and 
have been dimmed to correspond to host galaxy extinctions 
of Ar = 5.4 and 4.2 mag, respectively. For this we adopted a 
distance of 11 .7 Mpc for the host g alaxy of SN 1999em based 
on Cepheids (^Leon ard et al.|[2003h and 8.4 Mpc for the host 
galaxy of SN 2005cs based on the application of the Ex- 
panding Photosphere Method for SNe 200 5cs and 201 Idh in 
the same host galaxy (jVinko et al.ll2Q12l ). The total (host 
galaxy + Galactic) extinctions of Ay = 0.31 and 0.43 were 
adopt ed for S Ne 1999em and 2005cs, respectively, following 
Smartt et al.l ([2009) and a Galactic extinction of Ay = 0.53 
(|Schlegel et al.lll998l l towards M 82. Assuming an epoch of 



© 0000 RAS, MNRAS 000, 000-000 



4 Mattila et al. 



Table 2. Reported coordinates of the transients in M 82 discussed in this paper. All coordinates are J2000. 



Transient 


RA 


Dec 


Reference 


SN 2004am 


9:55:46.61 


+69:40:38.1 


Sineer et al. 


(2004) 


SN 2008iz 


9:55:51.55 


+69:40:45.792 


Brunthaler et al. 


(2009) 




9:55:51.55 


+69:40:45.788 


Beswick et al. 


^2009) 


MERLIN Transient 


9:55:52.5083 


+69:40:45.420 


Muxlow et al. 


(2009) 





Figure 1. A 10"xlO" region around SN 2004am in M 82 (North is up and East to the left). Kg-band image obtained with LIRIS on 
WHT on 2004 Nov. 25 (left), the subtraction between the LIRIS Kg-band images obtained on 2004 March 6 and 2004 Nov 25 (middle), 
and an ACS/WFC F814W image obtained on 2006 Mar 29 (right). The position of SN 2004am as determined from the subtracted LIRIS 
image is indicated in the HST image. 



explosion of 14 days before the SN was first observed by 
the LOSS was found to yield a good match with the post- 
plateau drop in the light curves of SNe 1999em and 2005cs. 
If the extinction toward SN 2004am is similar to that found 
for M 82-L by Lancon et al. (2008) then its absolute mag- 
nitude is fainter than that of the low- luminosity type II-P 
SN 2005cs. SN 2004am is only likely to be a normal type 
II-P SN if the extinction is higher than Ar 5. Also, the 
spectrum analysis indicates that SN 2004am has a low ejecta 
velocity and is similar to SNe 1997D and 2005cs. Further- 
more, SN 2004am is not visible in the late deep image from 
HST ACS at 241 days, allowing a conservative detection 
limit of / 20.3 to be set. The corresponding R-hand limit 
would be fainter than this due to the red R — I colour of SNe 
II-P at such epochs and the likely high extinction towards 
SN 2004am. Figure 3 shows this is too faint to be compatible 
with the nebular phase of SN 1999em, but is compatible with 
a SN that ej ected a low amount of ^^Ni such as SN 2005cs 
(|Pastorello et al. 2009). 

Comparison of the unfiltered SN photometry from 
Singer, Pugh & Li (2004) on 5 March with our ii^— band pho- 
tometry on 7 March yields a SN color of (mccD — i^) — +5.0. 
Typical R — K magnitudes around 100 days for II-P SNe are 
in the range 0.8 — 1.6 (Maguire et al. 2010), hence this sup- 
ports a host galaxy extinction of around Ay — 4 — 5. The 
NIR colours J - = 0.9 ± 0.2 and if - = 0.5 di 0.2 
also support a large extinction for SN 2004am. Given the 
evidence that low- luminosity SNe show a large NIR excess 
which changes rapidly at around 100 days (Pastorello et al. 
2009), the JHK colours are compatible with extinctions in 
the range Ay ^3 — 6. 



2.4 Super-star cluster M 82-L and the progenitor 
of SN 2004am 

The LIRIS i^s-band image obtained on 2004 Nov. 5 was 
aligned with the (10' x 10') HST/ACS W FC F814W (/- 
band) mosaic of M 82 (0.05 arcsec pixel" ^ ; iMutchler et alj 
2007). We measured the centroid positions for 13 point-like 
sources visible in both the images using the IRAF AP- 
PHOT package. A geometric transformation function in- 
cluding shifts, rotation and scaling in both x and y was de- 
rived for the pairs of coordinates using the IRAF GEOMAP 
task. The astrometric uncertainty was estimated from the 
RMS of the residuals from fitting the transformation func- 
tion to the data points. One source was discarded at this 
stage as an outlier. The uncertainties in x and y estimated 
this way were 0.17 and 0.18 arcseconds, respectively. Higher 
resolution images of the field around M 82-L and M 82-F are 
available from the HST/ACS HRC (0.025 arcsec pixel"^) 
and we have also employed the F814W filter image in mea- 
suring the clusters radial profile. These were taken on 2002 
June 07 and presented in Bastian et al. (2007). However the 
HRC images did not contain enough stars to astrometrically 
position SN 2004am, or align it accurately to our NIR im- 
ages. 

To measure the position of SN 2004am we first sub- 
tracted the Nov. 25 Xg-band image from the March 6 K^- 
band image using the ISIS2.2 package. The position of 
SN 2004am was then measured in the subtracted image 
making use of three different methods viz. centroiding, 2 
dimensional Gaussian fitting, and optimal filtering, within 
the IRAF APPHOT package. The mean of these coordi- 
nates was adopted as the best estimate for the position of 
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SN 2004am which was converted to the coordinates in the 
ACS WFC. The uncertainty in the position of SN 2004am 
determined as the standard deviation of the three differ- 
ent measurements was found to be neghgible compared to 
the RMS of the geometric transformation. We therefore find 
that the supernova is located 0.28=b0.25 arcsec North-East 
from the centroid position of the SSC in the ACS image 
(see Fig. 1). In the high resolution ACS HRC F814W image 
the SSC is resolved, but it has a fairly sharp core compared 
to t he nearby M 82-F a nd two other nearby, fainter clus- 
ters (Bastia n et al.|[2007l V The core has a half width at half 
maximum (HWHM) of approximately 0.03 arcsec, but the 
wings of the radial profile are broad and the half-width at 
zero intensity is 0.38 arcsecs. M 82-F by comparison has a 
broader core with a HWHM of 0.07 arcsec, but a similar 
spatial extent of approximately 0.4 arcsecs. While the best 
measured position for SN 2004am falls outside the half-light 
radius of the cluster M 82-L, the error box stills leaves it 
within the observed (and resolved) radial extent of the clus- 
ter flux. We conclude therefore that SN 2004am is spatially 
coincident with the SSC M 82-L and that its progenitor was 
most likely a cluster member. 



M 82-L has been extensively studied in the optical and 
the NIR , with the first attempt at an age determination of- 
fered bv lCallagher SmithI (|l999l ). Their spectra of M 82-L 
only included red wavelengths, resulting in an somewhat un- 
certain estimate of cluster age. However, they suggest an 
age of around 60 Myrs, similar to that obtained for the 
M 82-F, as these clusters display comparable Call triplet 
line strengths. Space Telescope Imaging Spectrograph (STIS) 
spectra of M 82-L were analysed by Smith et al. (20 061) with 
the spectral synthesis code starburst99 (|Leitherer et al.l 
I1999I I. A spectroscopic age for the cluster of 65I35 Myrs and 
extinction of E{B - V) = 1.87±0.39 was derived. The large 
error bars arise mostly from the fact the the [/-band was 
not covered and hence extinction and age become difficult 
to disentangle. 



The more recent study of Lancon et al. (2008') presented 
a medium resolution (R ^ 750) and high signal-to-noise 
0.8 — 2.4/im spectrum of M 82-L. They used the stellar pop- 
ulation synthesis code (PEGASE.2), with a library of ob- 
served red supergiant spectra to produce single stellar pop- 
ulation spectra with ages between 8 to 60 Myr. In all cases 
near-solar metallicity is assumed, along with a Salpeter ini- 
tial mass function (IMF). An impressive fit to the JHK 
spectra of M 82-L is obtained for a single stellar population 
of age 18lg^ Myr. The age of the cluster provides a strong 
constraint on the mass of the progenitor star, if we assume 
that the cluster formed coevally and t hat the cluster age is 
representative of that of the progenitor. ISmartt et aD (l2009ll 
employed the stellar evolutionary tracks of 'Eldridge ToutI 
(12004) to determine stellar masses of II-P progenitors. Us- 
ing these tracks (to allow a consistent comparison) we find 
that an age of 18lg^ Myr corresponds to lifetimes of stars 
of initial masses 12^ 3 M ^. The models in t he single stel- 
lar population analysis of Lancon et al. (2008) were those of 
[BressaneiaD (|l993) and the age-mass relationships in both 
codes are quite similar. 



3 SN 2008IZ AND THE MERLIN TRANSIENT 
- OBSERVATIONS AND RESULTS 

3.1 Observations and data reduction 

The sites of SN 2008iz and the MERLIN transient were 
observed on 2009 Jun 11 using the Near-Infrared Imager 
(NIRI) with the ALTAIR AO system on the Gemini-North 
Telescope, as part of our SN progenitor identification pro- 
gran|3- The f/14 camera was used, giving O.OS^^pixels over a 
51^^ X 51^^field of view. While typically in AO observations 
of SNe, the SN itself is used as a natural guide star (or a 
tip-tilt star together with a laser guide star), this is only 
possible for bright {V < 15 mag) SNe. Hence for these ob- 
servations we used the nearby cluster M 82-F as a natural 
guide star, which yielded a reasonable AO correction, de- 
spite M 82-F being an extended source. The image shown in 
Fig. 4 was reduced using standard techniques (sky subtrac- 
tion and masking of bad pixels) for NIR imaging with the 
IRAF. GEMINI package. 

As neither SN 2008iz nor the MERLIN transient were 
immedia tely visible at their positions from the radio obser- 
vations ([Eraser et al. |[2009h , we attempted to perform im- 
age subtraction to identify any faint transient source or low 
level variability at the site of either transient. We obtained 
pre-explosion observations from HST+NICMOS via MAST0 
to use as reference images for the image subtraction. For 
SN 2008iz, the images used were a 448s exposure with the 
NIC2 camera (pixel scale 0.075'7pix) in the F222M filter 
and an 1408s exposure in the F237M filter obtained using 
the same camera. For the MERLIN transient, the templates 
used were a 448s exposure in F222M and a 1408s exposure in 
F237M. All the NICMOS images were taken on 1998 April 
12; the pipeline- reduced products from the HST archive were 
used. 



3.2 Alignment and SN identification 

As the field of view of both the NIRI (51'' x 51'') and NIC- 
MOS (19" X 19") images are small compared to the angular 
size of M 82, we employed a bootstrap technique to obtain 
accurate astrometry for our data. We first identified SDSS 
sources in an F8141^-filter HST ACS mosaic (Mutchler et 
al. 2007) of M 82, which allowed us to calibrate the world 
coordinate system (WCS) of this image to a high degree of 
accuracy, and subsequently used this frame as an astromet- 
ric reference for all other images. 

The ACS mosaic of M 82 was downloaded from the 
HST MAST high level science products webpage^. This 
mosaic image has been corrected for the geometric distor- 
tion of ACS, and its internal astrometric accuracy is ex- 
tremely good, with residuals in the alignment of the input 
tiles to the mosaic of ^^0.1 pixels, corresponding to 5 mas. 
We identified 51 point sour ces which are in the Sloan cata- 
log ('Abazaiian et al. ]|2009h with an r magnitude < 21 mag. 
The coordinates of these sources were measured with the 
IRAF PHOT task. Using the pixel coordinates of the sources. 



1 GN-2009A-Q-32 

^ http://archive.stsci.edu 

^ http:/ /archive. stsci.edu/prepds/m82/ 
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Figure 2. Sections of ACS/HRC F814W images (North is Up and East to the Left) showing the SSCs M 82-L (to the south-west from 
the centre) and M 82-F obtained on 5 July 2004 (left) and 7 June 2002 (middle) and the subtraction between the two obtained using 
ISIS2.2 (right). The first two are shown with inverted intensity scale. Positive image subtraction residuals are visible in the subtracted 
image at the location of M 82-L but no evidence of a point source indicating that the SN has faded below the HST detection limit. 



together with the Sloan catalog celestial coordinates, we re- 
fit the WCS for the ACS mosaic using iraf CCMAP and 
CCSETWCS. A general transformation was allowed, consist- 
ing of translation, rotation, scaling and a polynomial term. 
Three sources were discarded at this stage as outliers. The 
rms error in fitting the new WCS was 48 and 41 mas in R.A. 
and Dec, respectively. 

Next, we identified 26 sources in common between our 
NIRI i^-band image and the ACS F814Ty-filter mosaic. We 
measured the pixel coordinates of these as before with iraf 
PHOT, and fitted the matched coordinate list with a general 
transformation using IRAF GEOMAP. We rejected two sources 
as clear outliers from the fit. The rms error in the fit was 
20 and 22 mas in R.A. and Dec, respectively. The NIRI 
image was transformed to match the ACS image, and the 
two frames compared by eye to verify the transformation. 

Using the Sloan- calibrated WCS, and using the values 
of R.A. and Dec. from radio observations given in Table[2]we 
identified the pixel coordinates of SN 2008iz and the MER- 
LIN transient in the ACS image. These coordinates were 
then transformed to the NIRI image using the transforma- 
tion determined previously together with the geoxytran 
task. 

Unfortunately, the site of SN 2008iz and the MERLIN 
transient are on different NICMOS images, and so each had 
to be aligned to the NIRI image separately. In each case, 
we matched between ten and twenty common sources, and 
transformed the NICMOS image to the pixel coordinates of 
the NIRI image using iraf geomap and geotran. 

For SN 2008iz, we used the isis2.2 package to convolve 
the F222M- and F237M-filter NICMOS images to match 
the PSF and flux level of the NIRI i^-band image. We then 
subtracted the convolved F222M and F237M images from 
the K image, and searched for residual flux at the site of the 
radio SN. As an additional test we subtracted the F222M 
image from the F237M image to see if a colour difference 
could produce a source. 

The results of these subtractions are shown in Figures 5 
and 6. As can be seen, there is a clear detection of a source 
co-incident with the radio coordinates of SN 2008iz in the 
subtractions between the NIRI and NICMOS images. The 
source in the difference image consists of positive flux, im- 



plying a brightening since April 1998. The source is clearest 
in the subtraction between NIRI K and NICMOS F222M 
images (which is unsurprising given the good match between 
the bandpasses of the filters, as shown in Figure[7|), although 
it is still significant in the F237 filter. No source is detected 
at the site of the SN in the subtraction between the F222M 
and F237M filters, indicating that there is not a point source 
with a large F222M - F237M colour at this location that 
could account for the source in the other subtractions. 

We repeated the same procedure for the MERLIN tran- 
sient, again finding a source at the radio position of the 
transient in the difference image. As the detection was not 
as clear as in the case of SN 2008iz, we repeated the convolu- 
tion and subtraction process using hotpant^, but obtained 
the same result. 

We measured the pixel coordinates of both sources in 
the difference images obtained with ISIS and the F222M 
template image using the optimal filter centering algorithm 
in IRAF PHOT. In the case of SN 2008iz, we find an oflFset 
of 0.68 pixels (34 mas) between the expected and measured 
position of the transient, while for the MERLIN transient 
there is a 1.62 pixel (81 mas) offset. While this is slightly out- 
side the combined rms error of the alignment (70 mas), we 
nonetheless consider the association convincing given that 
the transient source is detected at a relatively low signal to 
noise. We identified nine other point sources which appeared 
in a 12^^ x 12 ^^region in the diflFerence image for the MERLIN 
transient; from this we calculate a probability of a variable 
or transient source coincident to less than 81 mas as ^0.1 
per cent. 



3.3 Photometry of SN 2008iz and the MERLIN 
transient 

For both SN 2008iz and the MERLIN transient, we per- 
formed aperture photometry on the detected sources in the 
K-F222M difference images. PSF-fitting photometry was 
not used, as the PSF in adaptive optics imaging (i.e. the 
PSF to which the NICMOS image was convolved to match) 



http:/ /www. astro. washington.edu/users/becker/hotpants. html 
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Epoch (post explosion) 



Figure 3. Left-Top : A comparison of the SN 2004am zJ spectra obtained on 2004 Mar 6 (black) with the spectrum of SSC M 82-L 
from Langon et al. (2008) shown in red (along with a Gaussian smoothed spectrum overlaid in green). The pre-discovery spectrum of 
M 82-L was scaled and subtracted from the 2004 Mar 6 spectrum to leave a clean spectrum of SN 2004am. The major features are 
marked at the positions of the probable P-Cygni absorption troughs. Left-Lower : A comparison of the SN 2004am zJ spectrum with 
that of three other well studied II-P SNe near the end of the plateau phase. All spectra were dereddened and scaled to match the flux of 
SN 2004am at l.S/xm (see text). As the spectra were taken with a diversity of resolutions, all were broadened (with a Gaussian FWHM of 
40A) and rebinned (to 20Apix~-'^) to the lowest resolution for a meaningful comparison. The SN 2004am spectrum was dereddened with 
Ay = 3.7 and R = 2.4 as found for M 82-L by Langon et al. (2008). Right : The light curve of SN 2004am compared with the R-hsmd 
light curves of two type II-P SNe. The comparison light curves were scaled to a distance of 3.3 Mpc and dimmed to correspond to host 
galaxy extinctions of Aj^ = 5.4 and 4.2 for SN 1999em and 2005cs, respectively. The plateau points are unflltered CCD magnitudes from 
Singer, Pugh &: Li (2004), and the upper limit at 241 days is an estimated limit from the F814W HST image. An epoch of explosion of 
14 days before the flrst observations by LOSS was adopted for SN 2004am to yield a good match with the post-plateau drop in the light 
curves of SNe 1999em and 2005cs. 



varies as a function of position. To set the zero point for 
the photometry in the subtracted images, we used the pub- 
lished F222M magn itudes o f super-star clusters 6, 7, 8, 9 
and 11 from McCra dTet alJ |2003). We also attempted to 
use 2MASS sources in the NIRI field of view to set the zero- 
point, but unfortunately there was an insufficient number of 
isolated sources with reliable 2MASS i^-band magnitudes 
covered by the FOV. 

Using a photometric aperture with an 8 pixel radius 
for SN 2008iz gives a magnitude of F222M {^K) = 16.0 
for SN 2008iz. As a consistency check, photometry was also 
obtained from a subtracted image where the NIRI frame was 
scaled to match the flux units of the NICMOS template. 
Using the PHOTFLAM value in the FITS header of the 
latter to convert the counts measured for SN 2008iz to a flux, 

corresponding 



to i^^=16.14 in the CIT system (l Elias et al.lll982l ). We hence 
adopt i^=16.1 as the magnitude of SN 2008iz on 2009 June 
11. Varying the photometric aperture by a few pixels causes 
the photometry of SN 2008iz to vary by up to 0.2 mag. 
The measured magnitude is also dependent on the region 
used to measure the sky background at the SN location. The 
convolution used in the image subtraction will also affect the 
noise properties of the subtracted image. For these reasons, 
we adopt a conservative error of ±0.3 mag which is larger 
than would be implied by Poissonian statistics. 

For the MERLIN transient we used a smaller aperture 
with a 4 pixels radius to avoid including flux from an image 
artifact to the south of the transient in the difference im- 



gives 



. value of 1.34X 10"^'^ erg s"^cm"^J 



age (as can be seen in Fig. 6(c)). Using the zeropoint from 
photometry of nearby SSCs from McCradv et a l. (2003), we 
find a magnitude for the MERLIN transient of F222M = 
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Figure 4. A 50" x33" region of the K-band Gemini-N Altair/NIRI image of the M 82 nuclear regions (North is up and East to the left). 
The locations of SN 2008iz, the MERLIN transients and the SSCs M 82-L and M 82-F are indicated. 



17.5, with the same conservative error of 0.3 mag as for 
SN 2008iz. As a check on the effect of the small aperture on 
our photometry, we used the same 4 pixel radius to measure 
the magnitude of SN 2008iz, and found a value which differs 
only by 0.2 mag from that found through the 8 pixel aper- 
ture. As this difference is smaller than our errors, we do not 
regard this as a significant source of error. 

As a final test we used synphot to calculate the 
F222M — K colour of a range of black bodies between 100 
and 50 000 K. For temperatures hotter than 500 K, we find 
that the colour difference between the two filters is negligible 
(0.02 mag or less). While the difference does become large 
(>0.1 mag) for T<400 K, if the flux in SN 2008iz was coming 
from dust at this temperature, then it would be emitting in 
the mid-IR rather than the K band (^Mattila et al...2008a V 



4 DISCUSSION 
4.1 SN 2004am 

SN 2004am is the only SN ever discovered at optical/IR 
wavelengths in the prototypical starburst galaxy M 82. It 
occurred spatially coincident with the obscured SSC M 82-L 
and suffers from a host galaxy extinction of Kv — 5=bl which 
is consistent with the extinction derived for M 82-L. From 
the compiled photometry and our NIR spectra it is clearly 
a type II-P SN, which one would expect to come from a red 
supergiant progenitor star. It is likely a low luminosity type 
II-P SN and these have been proposed to come from stars 
toward the lower mass range that will produce core-collapse 
(Smartt et al. 2009). 




1.9 2 2.1 2.2 2.3 2.4 2.5 2.6 



Wavelength (microns) 
Figure 7. Transmission curves for HST (telescope + instrument 
+ filter) and Gemini (filter only) used for the pre- and post- 
explosion imaging of SN 2008iz and the MERLIN Transient. 



iLangon et al.l (|2008l l and iMcCradv GrahamI (|200/t ) 
show that the NIR spectra of M 82-L are dominated by 
the spectral features of red supergiant stars. The viral 
mass of the cluster i s estimated to be 4 di 0.6 x 10^ 
(McCr ady Graham| [2007) which agj rees with the mass es - 
timate from the cluster luminosity bv lLangon et al ] (|2008l l. 
Assuming the best fit age of around 18 Myr (SectlZ4]) for 
M 82- L, the single stellar population models of lLancon et al.l 
(2008) suggest there are likely to be of the order of 100- 
260 red supergiants in the cluster, depending on the initial 
mass function and lower mass limit used. Clearly there is a 
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explosion pre-explosion, transformed to using ISIS 

match NIRI image 




pre-explosion pre-explosion using ISIS 

Figure 5. Gemini post- and HST pre-explosion images of the site of SN 2008iz, together with difference images obtained with ISIS. In 
all cases the images are centred on the radio coordinates of the transient, with North up and East left. All images are 3" x 3'' . 



rich population of red supergiants as potential progenitors 
of type II-P SNe. 

In fact M 82-L is the most massive host cluster of a SN 
that has been accurately measured in the nearby Universe. 
Nearby SNe do not tend to be discovered coincident with 
the most massive unresolved clusters. In the volume lim- 
ited sample of 20 II-P SNe reviewed bv lSmartt et al.1 ([2009|) 
which have high-resolution pre-explosion images available, 
only two fall on compact cluste r s. Thi s may be a little sur- 
prising, however IPellerin et al.l (|2007l ) show that in NGC 
1313 75-90 per cent of the UV flux is produced by stars out- 
side the dense clusters. They propose that the late O and 
early B-type stars (that will eventually produce 8-30 Mq red 
supergiants) are diffusely spread througho ut a gala x y due 
to infant mortality of stellar clusters. Crock ett et al. 1^2008) 
have shown that the type Ic SN 2007gr exploded on the edge 
of a bright object in NGC 1058 and suggested that this may 



be a host cluster. This is one of 12 type Ib/c SNe studied by 
Eldridge et al. (in prep.) at high resolution and is the only 
one possibly associated with a stellar cluster, again support- 
ing an approximate figure of 10 per cent. 

The mass range that we suggest is appropriate for 
the progenitor (12^3 M©) in agreement with the masses 
of red super giants directly detected in nearby galaxies. 
Is mart t et al.l (200^) present a full list of a volume and time 
limited search for SN progenitors, and find the likely pro- 
genitor population of II-P SNe are red supergiants of initial 
masses 8-17 M0. The mass of the progenitor of SN 2004am 
is consistent with this mass range. The fact that it was in 
a cluster with such a rich red supergiant population fur- 
ther supports the idea that this stellar evolutionary phase 
directly gives rise to h ydrogen rich II-P SNe. However, we 
note that Era ser et al. I (^2011) suggest a lower mass range of 
between 8-9 Mq for the progenitors of sub-luminous Type 
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(a) Gemini+NIRI (Ks) post- (b) HST+NICMOS (F222M) (c) Subtraction of [6(b)] from [6(a)] (d) Subtraction Qf [6(b)] from [6(a)] 
explosion pre-explosion, transformed to using ISIS using hotpants 

match NIRI image 





(e) Gemini+NIRI (Ks) post- (f) HST+NICMOS (F237M) pre- (g) Subtraction of [6(f)] from [6(e)] (h) Subtraction of [6(f)] from [6(e)] 
explosion explosion, transformed to match using ISIS using HOTPANTS 

NIRI image 





(i) HST+NICMOS (F237M) pre- (j) HST+NICMOS (F222M) pre- (k) Subtraction of [6(j)[ from [6(1)] (1) Subtraction of [6(j)] from [6(1)] 
explosion, transformed to match explosion, transformed to match using ISIS using HOTPANTS 

NIRI image NIRI image 

Figure 6. Gemini post- and HST pre-explosion images of the site of the MERLIN transient, together with difference images obtained 
using ISIS and HOTPANTS. In all cases the images are centred on the radio coordinates of the transient, with North up and East left. All 
images are 3^^x ?>" . 



II-P SNe, and if SN 2004am was one of these events, it is 
also in agreement with the expected progenitor mass range 
within the la uncertainty. 



velocity of -20 000 km s"^ observed for SN 2008iz (Brun- 
thaler et al. 2010) was somewhat higher than the one for SN 
1993J f-15 000kms~^ Marcaide et al.[[ 19951 ) 



4.2 SN 2008iz 

In Fig. [SI we compare the radio luminosities of SN 2008iz and 
the two radio transients in M 82 (41.5+59.7 and 40.59+55.8) 
of unknown nature with the peak spectral radio luminosi- 
ties of SNe with different types. The luminosities of the two 
transients are comparable with those of type II-P SNe, such 
as SNe 1999em and 2004dj. The peak luminosity and time 
to reach the peak make SN 2008iz very similar to the well 
observed type lib SN 1993J thus suggesting that also their 
progenitor stars might be similar. However, the expansion 



At a distance of 3.3 Mpc, SN 2008iz has an absolute 
magnitude of K= -11.6+0.3 for zero host galaxy extinction. 
The first detection of SN 2008iz at radio wavelengths was 
in March 2008, with a best estimate for the explosion epoch 
of 2008 Feb. 18, implying that at the time of our NIRI ob- 
servation the SN was already —480 days old. Unfortunately 
there are few SNe with NIR light curves for comparison at 
this late phase. Ho wever, SN 1987A had an a bsolute magni- 
tude of K= -12.2 fSuntzeff & Bouchet" 1990) at +450 days, 
while at 300 d ays SN 2004et had an absolute magnitude 
of K = -13.35 (jMaguire et al.[[201Ql ). These magnitudes are 
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comparable to that of SN 2008iz with a modest amount of 
extinction (Ak < ^ 1 mag). 

iBrunthaler et all |201Q| 1 estimated the extinction to- 
wards SN 2008iz to be Ay = 24. 4 ma^, using the in tensity of 



the ^^CO ( J = 2 ^ 1) line from lWeiB et~all (l2QQlh. together 



with the CO to H2 ratio and the relation of Giiver Ozell 
(|2009l l between hydrogen column density and extinction. 
We find however, that whe n using the same data and re- 
lations as IBrunthaler et al.L we obtain a total extinction of 
Ay — 48.9 mag. This is exactly twice the value found by 
IBrunthaler et all , and we suggest that these authors may 
have neglected to multiply the value of N{H2) by two to con- 
vert to hydrogen nuclei before applyinjSi; the hydrogen nuclei 
to extinction relation of 'Giive r Ozell 

Ay — 48.9 mag implies a n ext inction in K of ^^5.5 mag 
(ICardelh. Clavton. M^islll989l l. Such a high extinction 
would put SN 2008iz at an absolute magnitude of K ^ -17 
mag, which is uncomfortably high for a SN at such a late 
phase. While a late-time near-infrared excess has been ob- 
served for several SNe, and attributed to dust formation in 
the cool, de nse shell located between the forward and reverse 
shocks (e.g. lFox et al.|[2009l ). a K-band excess of this magni- 
tude is quite unusual. Furthermore, such excesses are most 
commonly seen in Type Iln SNe, which are intrinsically rare. 
For example, adopting a distance of 18.6 Mpc (assuming the 
Virgo-hGA+Shapley corrected recession velocity from NED 
and Ho = 70 km s"^ Mpc" ^) for the Type Iln SN 1998S its 
absolute X-band magnitude (|Pozzo et al.ll20o3 ) at an epoch 
of ^480 days was ^ -17.5. A more likely explanation of the 
apparent discrepancy between the extinction for SN 2008iz 
as derived from the H2 column density, and that expected 
for a SN at that phase is that much of the H2 lies behind 
the site of SN 2008iz. If the SN was located behind one fifth 
of the total H2 column, then its derived extinction would be 
closer to Ay ^ 10 mag, bringing the late time magnitude 
into agreement with that of SNe 1987A and 2004et. 



4.3 MERLIN transient 

One of the p ossible explanat i ons ad voca ted for the MERLIN 
Transient by IMuxIow et al.l (|2010l l and Ijoseph et al ] (|201lh 
was that it was an extra-galactic microquasar. Microquasars 
are the stellar analogs to quasars, and are believed to arise 
from a compact object accreting matter from a stellar com- 
panion, and in the process forming a relativist ic jet. They 
typically display strong and variable radio and X-ray emis- 
sion. In Figure [8] we compare the luminosity of the MERLIN 
transient with the range of luminosities that the Galactic 
microquasar Cyg nus X-3 has been ob served to reach in ma- 
jor fiare states (Waltman et al."l995'). Only if Cygnus X-3 
would increase its highest observed radio luminosity by a fac- 
tor of ^10, the MERLIN transient would be compara ble in 
brightness. We note that recent Iv .Batei at et al.l (|2012l l have 
proposed that three rapidly variable radio sources within 
the nuclear regions of the ultraluminous infrared galaxy 
Arp 220 having their radio luminosities a factor ^100 times 
larger than the luminosity of the Merlin transient could be 
associated with highly beamed microquasars. The appar- 
ent magnitude of the MERLIN transient is K—17.5, cor- 
responding to an absolute magnitude of K= -10.2 if as- 
suming no host galaxy extinction. Such a bright absolute 
magnitude would make the MERLIN transient the bright- 



10^ 



10 



10 



bD 10 



10 




10 



CygX-3 



O lb 

□ Ic 

A 11 

41.5+59.7 

40.59+55.8 

+ SN2008iz 
MERLIN transient 



10" 



I 1 1 
10^ 



10 



1 1 1 
10^ 



tp^^ [days]x[v/5GHz] 



Figure 8. Peak spectral radio luminosity versus time to reach 
the peak (normalized to 5 GHz) for different types of CCSNe. The 
four M 82 transients have been added to allow comparison with 
historic radio SNe. SN 2008iz is represented by a cross, whereas 
the rest of the transients, for which the peak epoch is unknown, 
are represented as dash-dotted (41.5+59.7), dotted (40.59+55.8) 
and dashed (MERLIN transient) lines. The gray region repre- 
sents the range that the galactic microquasar Cygnus X-3 has 
been observed to reac h in major flare states at radio wavelengths 
jWaltman et al.lll99dV 



est known example of a microquasar by quite some margin; 
for comparison, the Galactic micro quasar SS433 has an ab- 
solu te magnitu de which varies between K=-5.7 and 
6.4 (iKodaira. Na kada Backman I QSSi : iBlundeh Bowled 
I200I 

The MERLIN transient was first detected in the ra- 
dio 40 days prior to the e poch of the Gemin i ima ge. From 
the templates presented in Mattila & Meikle (2001) at +40 
days, we would expect a normal, unextinguished SN in M 82 
to have an apparent magnitude of K^9.6. Glearly the ob- 
served magnitude of the MERLIN transient appears difficult 
to reconcile with a SN. If, however, the first detection of the 
MERLIN transient did not mark the explosion epoch, but 
was instead late time radio emission, then the faint abso- 
lute magnitude is less puzzling. Similarly, ^70 magnitudes 
of extinction in V would bring the observed and expected 
K magnitudes into agreement. However, even with high ex- 
tinction or an earlier explosion epoch, the absence of any 
observed e xpansion at radio wav elengths (as was seen for 
SN 2008iz, IBrunthaler et al.ll2010l l poses a significant prob- 
lem. 

The absolute magnitude of the MERLIN transient is 
comparable to the quiescent magnitude of Luminous Blue 
Variables (LBV). LBV eruptions, also termed "supernova 
impostors" can reach an absolute magnitude of ^ -12 to 
-16 in V. The SN impostors SN 2009ip and UGG 2773 
OT2009-1 have upper limits on their 8 GHz radio luminos- 
ity of < 1.3 X lO^^erg s"^ Hz"^ and < 2.6 x 10^^ erg s"^ 
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Hz ^ respectively ('Foley et al.l l201lh . These limits are an 
order of magnitude higher than the peak radio luminosity 
of the MERLIN transient, and hence the luminosity of the 
latter is not inconsistent with the limited information on the 
properties of SN impostors at radio wavelengths. 

The NIR magnitude of the MERLIN transient is also 
reminiscent to that of a nova at peak. Most common no- 
vae fade rapidly af ter outburst, but fr om the template light 
curves presented in lStope et al either an F or D-class 

nova could still be bright at ~-|-40 days. However, classical 
novae are not usually strong radio sources, and the peak 
of the radio light curve tends to be later than the optical 
maximum, hence the MERLIN detection would appear in- 
congruous with this scenario. 

The nature of the MERLIN transient remains elusive, 
and on the basis of the limited data available, we regard 
an extremely bright extragalactic microquasar as the most 
plausible scenari o, perhaps from a high-mass X-ray binary 
such as LS 5039 (jClark et al.ll200lh . but with a high ratio of 
radio to X-ray flux. A bright extragalactic mic roquasar wa s 
propo sed a s the most likely ex planation by IMuxIow et al.l 
(|201Q| ) and I Joseph et al.1 (|201lh . We concur with that con- 
clusion although it would mean that the NIR luminosity is 
a factor of about 30 higher than Galactic microquasars. 



4.4 Comparison with the expected SN rate 

SNe 2004am and 2008iz are still the only supernovae de- 
tected at optical/NIR wavelengths in M 82. This is per- 
haps somewhat surprisin g ^ven the SN rat e estimates of 
aroun d 7-9 per century (i Fenech et al. 1 l2008l : iFenech et al. 
I2OIOI), the small distance and the appeal of this galaxy 
as a target in SN searches. It is very likely that extinc- 
tion has hampered the discovery of recent SNe in the past 
decade (for estimates of extinctions towards the radio SNRs 
of M 82 see iMattila k Meikld [2OO1I ) . Over the last -20-30 
years during which M 82 has been regularly monitored at 
radio wavelengths, the detection of two confirmed SNe (SNe 
2004am and 2008iz) and the two radio transients (41.5+59.7 
and 40.59+55.8) with a possible SN origin is in reasonable 
agreement with the expectation. However, the discoveries 
of SNe 2004am and 2008iz have important lessons for at- 
tempts to find heavily extinguished SNe in more distant 
starburst galaxies. While SN 2008iz was reco vered at a com- 
parable magnitude to the two SNe found by iKankare et al.l 
(120 12j) in th eir NIR LIRG SN search, it is of note that 
iFraser et al. ] (Hooa) originally reported a non-detection for 
both SN 2008iz and the MERLIN transient. We have only 
recovered these transients with a careful analysis using a 
posteriori knowledge of the source positions from radio ob- 
servations. Without these data, it is likely that SN 2008iz 
would have re mained undiscovered in the NIR. 

Recently iMattila et all (|2012l ) studied the fraction of 
CCSNe missed by rest-frame optical SN searches and found 
an average local value of ^-^20% increasing to ^-^40% of CC- 
SNe missed at z ^-^1-2. These estimates highlight the need 
for a better control of the SN activity in the obscured en- 
vironments such as the nuclear regions of nearby starburst 
galaxies and LIRGs. Unless properly accounted for such sys- 
tematic effects can domina te the uncertain t ies in the CCSN 
rates at high redshift fe.g.. iMelinder et"al ] (|2012l )). 



5 CONCLUSIONS 

SN 2004am was the first optically detected SN in M 82 when 
discovered by LOSS (Singer et al. 2004), and remains the 
only transient discovered in the optical. We show that it is 
most likely to have been a sub-luminous, highly reddened 
Type II-P event coincident with the obscured nuclear SSC 
M 82-L. From the cluster age we infered a progenitor mass of 
I2I3 M0which is within the uncertainties in agreement with 
the expected progenitor mass range for such events. Making 
use of high spatial-resolution K-band imaging we detected 
NIR counterparts for both SN 2008iz and the MERLIN tran- 
sient previously detected only at radio wavelengths. Our 
late-time K-hand magnitude rules-out an extremely high 
extinction towards SN 2008iz. The nature of the MERLIN 
transient still remains elusive, an extremely bright micro- 
quasar in M 82 being the most plausible scenario. 
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